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ABSTRACT 

We present and test TesseRACt, a non-parametric technique for recovering the concentration of sim¬ 
ulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos 
that are axisymmetric, triaxial, and contain substructure and compared to traditional least-squares 
fitting as well as two non-parametric techniques that assume spherical symmetry. TesseRACt recovers 
halo concentrations within 0.3% of the true value regardless of whether the halo is spherical, axisym¬ 
metric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry 
can return concentrations that are systematically off by as much as 10% from the true value for non- 
spherical halos. TesseRACt also performs significantly better when there is substructure present 
outside 0.5/^200* Given that cosmological halos are rarely spherical and often contain substructure, 
we discuss implications for studies of halo concentration in cosmological N-body simulations including 
how choice of technique for measuring concentration might bias scaling relations. 

Subject headings: galaxies: fundamental parameters — galaxies: halos — methods: numerical 


1. INTRODUCTION 

There has long been a disconnect between the way we 
describe a dark matter halo and the reality of that dark 
matter structure. In theoretical terms, we think of a halo 
as a smooth, isotropic, virialized, usually spherical, typi¬ 
cally uniformly rotating, distribution of mass that obeys 
a distinct radial distribution. However, while halos in 
cosmological simulations and observations may conform 
to these in a statistical sense, any one halo is not really 
any of these things. Halos are triaxi al, anisotropic , and 
contain signific ant substructure (e. g. Jing & Suto 2002 


Gao et al.|2004||Despali et al.|20i4| Groener V Goldberg 

2014|). Despite this well known fact, many analysis tech¬ 

niques try to extract halo properties by imposing one or 
more of these assumptions. For example, any proced ure 
which fits a radial profile, be it Hernqnist (H ernqnist 
19901), NFW ([Navarro et al.|[l996), or Einasto (|Einasto 
&; Haud| 1989), to a halos mass distribution assumes the 
halo is both spherical and smooth. In an era of simu¬ 
lations quickly surpassing 10 9 particles, there is a need 
for physically motivated analysis techniques that do not 
impose constraints on what a halo should look like. 

Halo concentration is a particularly useful statistic for 
characterizing halos. Since halos that gain the major¬ 
ity of their mass at earlier times (when the mean den¬ 
sity of the universe was higher) should be more compact, 
concentration is believed to encode a great deal of infor¬ 
mation about halo formation and growth. There have 
been numerous studies on the relationships between halo 
concentration in cosmological simulati ons and h alo mass 
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However, claims are often conflicting and the majority o 


However, claims are often conflicting and the - majority of 
techniques used to measure concentration fall victim to 
the above assumptions. 

We propose a non-parametric method for estimating 
halo concentration using Voronoi tessellation that we dub 
Tessellation based Recovery of Amorphous halo Con¬ 
centrations (TesseRACt). Section ^[briefly describes 
Voronoi tessellation and outlines TesseRACt, ^summa¬ 
rizes several tests, and Q summarizes our findings and 
describes studies that can benefit from TesseRACt. 


2. THEORY/BACKGROUND 
2.1. Measuring Concentration 

The concentration parameter is traditionally defined in 
an NFW halo as 


Cnfw 


R200 
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where R200 is the radius enclosing a mean density that 
is 200 times the critical density of the universe and R s 
is the scale radius. Since R200 can be easily found, con¬ 
centration is typically obtained by fitting Eqn. [2] to the 
radially enclosed mass profile to find R s . 
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Although this is simple in theory, fitting even spheri¬ 
cal halos without substructure can be difficult. Fits can 
be highly sensitive to resolution at the center, deviation 
from the expected NFW po wer laws, and choice of bin¬ 
ning (See|Prada et al.|2012|). These effects can be allevi¬ 
ated in practice by avoiding fitting altogether. Instead, 
the unknown profile parameters are related to other halo 
properties that can be robustly measured. For example, 
if Rhaif (th e radius enclosing half the mass) and R200 of 
a halo are known, it is possible to numerically solve 


1 _ In [(Rs + Rhalf)/Rs\ ~ Rhali/(Rs + Rhalf) 

2 In [(Rs + R20o)/Rs] — R200/(Rs + ft^oo) 
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for R s . 

This can be done for any two independent halo prop- 
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from NFW and yield more accurate results than fitting 
for dense halos with under sampled central regions, even 
these techniques still assume that halos are spherical and 
do not contain substructure. 

In principle, more accurate concentrations for non- 
spherical halos could be obtained by fitting to triaxial or 
ellipsoidal bins. Such techniques have been found to pro¬ 
vide more accurate mass estimates in both simulations 
and observations, but generally assume that the axis ra- 
tio and alignment remains constant throughout the halo 
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comi ng less and less spherical as yon look deeper in the 


halo (Allgood et al. 2006 Vera-Ciro et al. 2011). This 
makes measurements assuming a constant shape depen¬ 
dent upon the location within the halo at which shape 
is measured. In addition, despite allowing for more free¬ 
dom in halo shape, non-spherical fitting is still victim to 
the same caveats as spherical fitting and relies on the 
additional measurement of halo shape. 

While both non-spheroidal binning and non- 
parametric spherical techniques have advantages, 
neither is completely free of assumptions or can handle 
substructure. However, using Voronoi tessellation, we 
can construct a technique that does not rely on fitting, 
does not make any assumptions of spherical symmetry, 
and allows for substructure. 


2.2. Tessellation Based Concentration 

Given a set of seed points {pi,... ,p n } in some space, 
Voronoi tessellation divides the space between the seeds 
such that each seed is the closest seed to its Voronoi 
region. In this way, each seed (pi) has a correspond¬ 
ing Voronoi region of volume W encompassing all points 
in space which are closest to that seed. Voronoi tes¬ 
sellation has been used to non-para metrically identify 
galax y clusters in galaxy surveys (Soares-Sa ntos et al . 
201 1|), identi fy dark m atter halos (Neyrinck et al.||2004p 
and voids (Neyrinck 2008) in cosmological simulations, 
and imp rove the treatme nt of hydrodyn amics in simu¬ 
lations (Mocz et al. 2013| |Hopkins 2014). We take this 
one step further. Once halos are identified in cosmo¬ 
logical simulations (either by friends-of-friends, spherical 
over-density, or tessellation), the additional information 
provided by the particles’ associated volumes can be used 
to derive halo properties (like concentration) without im¬ 
posing any additional functional form. 

To determine concentration from the particle volumes, 
a profile is constructed that describes how mass scales 
with volume rather than radius. For a particle Pi with 
mass mi and a Voronoi volume V, the volume V enc ^ ‘en¬ 
closed’ by that particle is taken to be the sum of all 
particle volumes which are smaller than Vi or 
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Similarly, the mass M enc ^ ‘enclosed’ by a particle pi is 
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3=0 
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Each particle can then be assigned a theoretical ‘radius’ 
R[ = (3V„c,i/47r) 1 /3 that is defined as the radius the 
particle would be at if V enc ^ were spherical. The result 
is a volume based mass profile M enc (R'). Naively, the 
volume based concentration could then be defined as 


Cvol 



(6) 


where Woo is the densest volume containing an average 
density that is 200 times the critical density of the uni¬ 
verse and V s is some scale volume. If the theoretical radii 
(R '200 and K) associated with these volumes converged 
to the corresponding physical radii (R 200 and R s ) in the 
case of a spherical halo, c vo i would equal c n f w and this 
would be sufficient. However, this is not strictly true. 

For even a spherical halo, the relationship between a 
particle’s physical radius Ri and theoretical radius R[ is 
not 1:1. Due to the intrinsic scatter in the inter-particle 
spacing at a given physical radius, particles with slightly 
larger/smaller volumes will be scattered to larger/smaller 
theoretical radii. Because the density of particles is al¬ 
ways greater toward smaller physical radii, it is more 
likely for particles inside Ri to be scattered to larger 
theoretical radii. As a result, there will then be system¬ 
atically fewer particles considered ‘enclosed’ by particle 
Pi and R[ will be systematically lower. 

In order to correct for this and preserve the same nu¬ 
merical values for c vo i and c n f w in the case of spherical 
halos, the volume based concentration is defined as 


c V oi — (d 



( 7 ) 


/3 = 0.8062 and a = 1.0417 were obtained by fitting 
to measurements of Woo and V s for 10 spherical halos 
with known concentrations between 5 and 7 0. W hile this 
treatment is simplistic, tests performed in §3.1| indicate 
that it should be sufficient for most studies. 

Eqn. [7] can then be simplified in terms of the radii a 
halo would have if it were spherical. 


<8) 

Note that since Woo is defined as the densest volume con¬ 
taining an average density that is 200 times the critical 
density of the universe, R' 200 can be found directly from 
the Voronoi volumes. However, V s and R' s are somewhat 
arbitrary in the absence of fitting. Instead, we choose to 
define the scale radius and volume in terms of a quan¬ 
tity we can measure. We adopt Waif = 47ri2{f alf /3, the 
densest volume enclosing half the mass. 

Once R' haU is known, the corresponding R' s for a spher¬ 
ical halo can be found by numerically solving Eqn. [3j 
This relationship yields unique concentrations for each 
.Rhaif when 1 < c < 100. Above c = 100, Rhaif /^200 
becomes somewhat degenerate about 0.1. 

The TesseRACt procedure for finding c vo i is then: 
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1. Run Voronoi tessellation to determine the volumes 
and densities at each particle. 

2. Rank particles in order of decreasing density (in¬ 
creasing volume). 

3. Calculate enclosed mass (M enc ) and volume (V enc ) 
at each particle by summing the masses and vol¬ 
umes of ‘denser’ particles. 

4. Calculate mean enclosed density at each particle 
((Penc) = M enc /V enc ') 

5. Find V 200 and M 200 (the volume and mass enclosed 
by the particle at which the mean enclosed density 
reaches 200 times the critical density of the uni¬ 
verse) 

6. Find Vhaif (the volume enclosed by the particle at 
which the enclosed mass reaches M 200 / 2 ) 

7. Calculate R 200 and R'h a u f rom R 200 and Vhaif 

8. Numerically solve Eqn. [ 3 ] for R' s 

9. Calculate c vo i from Eqn. [8] 


3. TESTS 


For the tests that follow, we compare our volume- 
based non-parametric technique TesseRACt to both 
least-squares NFW fitting, as well as two non-parametric 
techniques that assume spherical symmetry. The first in¬ 
fers the scale radius from the Rhaif and R 200 in the same 
fashion as Eqn. [3J but using the physical particle radii. 
The second uses the relationship between the pe ak cir¬ 
cular velocit y and circular velocity at R 200 from |Prada 
et al.l (120121), 


^200 


0.216c 


1 1/2 


ln(l + c) - c/( 1 + c)_ 


(9) 


For each test, we run all four techniques on isolated N- 
body halos generated by sampling spherical NFW pro¬ 
files of known concentration. Unless otherwise stated, 
each halo has 7V part = 1 x 10 6 particles, M 2 00 = 1 x 
10 12 M© , and R 2 00 = 200 kpc. The halos only differ in 
po and R s . 

To ensure that insufficient sampling of the profile does 
not affect the accuracy of the least squares fitting, fits 
are bounded on the lower end at O.O 5 R 200 and R 200 on 
the upper end. 

Errors on concentration measurements were deter¬ 
mined by running each technique on 10 different real¬ 
izations of each test halo generated by using a different 
random number seed. In the tests below, the mean con¬ 
centration returned by each technique is plotted with the 
standard deviation across the different realizations. 


3.1. Concentration 

We first tested the performance of each method on ha¬ 
los with concentrations of 5, 10, 25, and 50. Results are 
shown in Figure [lj 

As these are idealized halos, all of the tested techniques 
return accurate measures of concentration as expected. 
For the idealized halos, TesseRACt is the most accurate 
for all but the c = 5 halo for which the measured con¬ 
centration is still within 2% of the correct value. Of the 
techniques which assume spherical symmetry, all three 
perform similarly with the maximum circular velocity 



Fig. 1.— Top: Concentration measured by TesseRACt and dif¬ 
ferent spherical techniques as a function of known halo concen¬ 
tration. Bottom: % relative residuals, 100*(measured-true)/true. 
The black line is the known concentration, blue is TesseRACt, red 
is least-squares fitting to Eqn. [ 2 ] green is the half mass technique 
from Eqn. [3] assuming spherical symmet ry, and magenta i s the 
maximum circular velocity technique from Prada et al. 1(2012). All 
of the techniques recover the correct value to within 2 %, regardless 
of concentration. 


and fitting techniques returning marginally more accu¬ 
rate concentrations at the low and high ends respectively. 
The accuracy of the techniques which assume spherical 
symmetry does not appear to be overly concentration 
dependent. However, as the modified concentration def¬ 
inition used for TesseRACt in Eqn. [8] was fit to mea¬ 
surements across a range of concentrations, TesseRACt 
is slightly more accurate for halos with intermediate con¬ 
centrations. All of the techniques tested were consistent 
across the 10 halo realizations with standard deviations 
of < 0.6%. 


3.2. Halo Shape 
3.2.1. Oblateness 

To test how each method performed on oblate halos, 
the c = 10 halo from above was squeezed along the 2 axis 
using the volume preserving transformation 
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for 7 different values of ellipticity where e Q biate = b/a and 
a = 1. Results are shown in Figure [2j 

For all but the least spherical halos, TesseRACt consis¬ 
tently recovers concentrations within 0.5% of the correct 
value, with a minor dependence on halo shape (overes¬ 
timated by 0.4% for e Q biate = 0-3 vs. underestimate by 
0.02% for spherical halos). However, the performance 
of the techniques which assume spherical symmetry is 
highly dependent on halo shape. Fitting performs the 
worst overall, with increasingly underestimated values 
for the least spherical halos (10% for e Q biate = 0-3). The 
performance of the two non-parametric spherical tech¬ 
niques is dependent on halo shape in a more complicated 
manner. 

For decreasing values of e 0 bi a te (decreasing spher¬ 
ical symmetry), both non-parametric spherical tech¬ 
niques have moderately increasing positive residuals until 
around e Q biate = 0.6 to 0.7 where the residuals begin to 
rapidly decrease resulting in underestimates of 6 and 8% 
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Oblate Ellipticity 



Fig. 2. — Performance of TesseRACt vs. spherical techniques 
as a function of oblate halo ellipticity. While TesseRACt recovers 
concentration with 0.5% of the correct value for all halo shapes, 
techniques assuming spherical symmetry do not and underestimate 
the concentration of the least spherical halos by up to 10%. 


for the half-mass and peak-velocity techniques respec¬ 
tively. These dependencies arise because the deformation 
of the halo flattens the radial mass profile. As fitting uses 
the entire profile, the concentration is underestimated for 
all non-spherical halos. However, non-parametric tech¬ 
niques use only two points in the profile (R200 and some 
inner radius) and their performance will depend up how 
these parts of the profile are affected by the transforma¬ 
tion. R200 continually increases as compression of the 
halo edge decreases the density at the edge of the profile. 
The inner radius will increase as well, but at a slower 
rate resulting in an overestimation of the concentration 
until the edge of the halo is compressed to a size com¬ 
parable to the inner radius. Once this occurs, the inner 
radius increases more rapidly than R200 resulting in un¬ 
derestimates. The precise value of e Q biate at which this 
occurs will depend on the radius that a particular non- 
parametric technique utilizes. Since the velocity peaks at 
a smaller radius than the half-mass radius for an NFW 
profile, the technique which uses the peak velocity is less 
sensitive to this effect. 

The test halos here have the same ellipticity at all 
radii by design. However, real halos are found to b e de- 
creasingly spherical at their centers (Allgo od" et al.| 2006 


Vera-Ciro et al. 2011). While TesseRACt would not be 
affected, it Is likely that the spherical techniques would 
produce less accurate concentrations. This would be par¬ 
ticularly pronounced for fitting due to sensitivity to the 
inner profile slope. 


3.2.2. Prolateness 

To test how each method performed on prolate halos, 
the c = 10 halo from above was squeezed along it’s z and 
y axes using the volume preserving transformation 
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for 7 different values of ellipticity where e pro iat e = b/a = 
c/a and a = 1. Results are shown in Figure p] 

The performance across all techniques is similar to the 
above tests for oblate halos. However, there does ap¬ 
pear to be a stronger dependence on ellipticity at inter¬ 
mediate values for prolate halos than for oblate halos for 


Fig. 3. — Performance of TesseRACt vs. spherical techniques as 
a function of prolate halo ellipticity. TesseRACt recovers concen¬ 
trations to within 0.3% regardless of halo shape. The performance 
of techniques assuming spherical symmetry are highly dependent 
on halo shape and return less accurate halo concentrations by over 
10%. Overall, the spherical techniques are less accurate at recov¬ 
ering concentrations for prolate halos than oblate halos. 


techniques which assume spherical symmetry. This is be¬ 
cause the prolate transformation flattens the radial mass 
profile to a greater degree than the oblate transforma¬ 
tion. 


3.2.3. Triaxiality 

We then tested the performance of each method on 
triaxial halos (See Figure |4|. For 9 di fferent values o f 

triaxiallity (T = b,F h a 2 , for c < b < a, Franx et al. 1991), 
the c = 10 halo from above was squeezed along the z and 
y axes using the volume preserving transformation 
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where a = 1, c/a = 0.3, and b = yfo 2 — T(c 2 — a 2 ). 



Triaxiality 

Fig. 4. — Performance of TesseRACt vs. spherical techniques as 
a function of halo triaxiality. The dotted colored lines denote the 
fiducial value returned by each technique in the case where T — 1 
(an oblate halo with e 0 bi a te — 0-3) and the residuals plotted here 
are with respect to this value. TesseRACt recovers concentrations 
within 0.2% of the fiducial value for all triaxialities. The techniques 
which assume spherical symmetry display moderate dependence on 
triaxiality, but do not return concentrations that deviate from the 
fiducial value by more than 4%. 


TesseRACt out performs all spherical techniques at 
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< 0.2% for all triaxialities. The spherical techniques re¬ 
turn concentrations that are only slightly less accurate 
overall (<4%) than in the case where T = 1. While 
TesseRACt’s performance does not depend on triaxial- 
ity, the accuracy of the techniques which assume spher¬ 
ical symmetry does to a small degree. The two non- 
parametric techniques become less accurate for lower val¬ 
ues of T, while fitting only exhibits this behavior until 
T ~ 0.3, where this trend is reversed. However, since the 
residuals never exceed 4%, it is reasonable to conclude 
that only the ratio between the largest and smallest halo 
axes appears to have a significant impact on the accuracy 
of the concentration estimate. 

3.3. Halos with Substructure 

To test the ability of TesseRACt to recover concentra¬ 
tions in the presence of substructure, idealized subhalos 
of varying mass and concentration were added to the 
c = 10 halo at varying radii. In each case the mass of 
the subhalo was set by downsampling the test halo from 
the previous section with the desired concentration. The 
size of the subhalo was then scaled such that the mean 
density within 0.2R 2 oo was half that of the parent halo. 

3.3.1. Substructure Mass 

Figure [5] shows the results from varying the mass of the 
subhalo across 1-20% of the parent halo’s M 2 oo- The sub¬ 
halos had a concentration of 50 and were placed 0.5R 2 oo 
from the center of the parent halo. Since subhalos were 
scaled to have a constant central density, increasing the 
mass of the subhalo also increases its size. 



Fig. 5. — Performance of TesseRACt vs. spherical techniques 
as a function of subhalo mass. The concentrations returned by 
all techniques become less accurate as the mass of the subhalo 
increases. However, at the most massive end, TesseRACt is three 
times as accurate as techniques assuming spherical symmetry. 


As expected, all techniques become less accurate as 
the subhalo increases in mass and size, but TesseR¬ 
ACt performs significantly better. For the most mas¬ 
sive subhalo, 20% the mass of the parent, TesseRACt 
overestimates the concentration by 8%, while techniques 
assuming spherical symmetry underestimate concentra¬ 
tions by up to 30% of the true value. TesseRACt overes¬ 
timates concentrations because particles in the subhalo 
have small enough volumes that they are assumed to be 
within the inner parts of the parent halo and thus result 


in a smaller half-mass volume. The techniques assuming 
spherical symmetry underestimate concentrations in the 
presence of substructure because the subhalo contributes 
mass to the outer profile. 

3.3.2. Substructure Radius 

Figure [6] shows the results from varying the distance of 
the subhalo from the center of the parent halo from 0.05 
to 0.5R 2 oo- The subhalos each had a concentration of 50 
and 1% the mass of the parent halo. 



Substructure Position 

Fig. 6. — Performance of TesseRACt vs. spherical techniques as 
a function of subhalo location in the parent halo. All techniques 
perform poorly when the subhalo is close to the center of the parent 
halo. However, while TesseRACt becomes more accurate as the 
subhalo is placed further out, spherical techniques do not. 


At small radii, the subhalo almost coincides with the 
center of the parent halo, causing all techniques to re¬ 
turn higher concentrations than they would had there 
been no substructure present. As the subhalo is placed at 
larger radii, the concentrations returned by all techniques 
become lower. However, while concentrations returned 
by techniques assuming spherical symmetry continue to 
drop past the true value as additional mass is added to 
the outer profile resulting in underestimates by up to 
20%, TesseRACt becomes more accurate, underestimat¬ 
ing the concentration by only 0.7% at 0.75R 2 oo* This oc¬ 
curs because, when the subhalo is placed at larger radii 
where the density of the parent is lower, the densities 
assigned to the subhalo particles using tessellation be¬ 
come lower and contribute less to the half-mass volume. 
The radius at which the spherical techniques change from 
overestimating to underestimating depends upon which 
part of the mass profile is used to calculate concentration. 

3.3.3. Substructure Concentration 

Figure [7] shows the results from varying the concen¬ 
tration ofthe subhalo from 5 to 50. The subhalos were 
placed 0.5R 2 oo from the center of the parent halo and 
had 1% the mass of the parent. 

For subhalos of equal or greater concentration than the 
parent (c su b > 10), TesseRACt is more accurate than the 
spherical techniques (< 3% vs. <20%). While the accu¬ 
racy of the spherical techniques shows little dependence 
on concentration in this regime, TesseRACt’s accuracy 
changes with concentration. Since TesseRACt places no 
constraints on where the center of the halo is, the mea- 
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Fig. 7. — Performance of TesseRACt vs. spherical techniques as 
a function of subhalo concentration. TesseRACt is more accurate 
than the spherical techniques overall, but exhibits some depen¬ 
dence on the concentration of the subhalo. The spherical tech¬ 
niques only exhibit a slight dependence for the subhalo that is less 
concentrated than the parent, where the non-parametric spherical 
techniques become more accurate. 
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Number of Particles 

Fig. 8. — Performance of TesseRACt vs. spherical techniques 
as a function of particle number. Dotted colored lines denote the 
fiducial value returned by each technique in the case of the full 
resolution, prolate halo. Residuals were computed with respect to 
this value. TesseRACt , fitting, and the half-mass technique are 
much less dependent on particle number than the technique which 
uses the peak velocity. 


sured concentration is essentially an average of the con¬ 
centrations of the two halos present, the parent and the 
subhalo. As a result, the s ame dependence on halo con¬ 
centration that was seen in §3.1| is also present for subhalo 
concentration. 

The case of a subhalo that is less concentrated than 
the parent halo (c su b = 5) is less physical. In order to 
maintain a central density half that of the c = 10 par¬ 
ent halo with 1% of the mass, the c = 5 halo must be 
scaled to less than half of its original size. This means 
that the substructure contributes to a narrower region 
of the radial mass profile, resulting in lower estimates 
for R 200 and M 200 that are closer to those for the par¬ 
ent halo alone. For fitting, which uses the entire profile, 
the result is a significantly underestimated concentration 
(20%). However, the non-parametric techniques, which 
are not dependent on the whole profile, are significantly 
more accurate in this case. The half-mass technique is 
slightly more accurate (2%) than the peak velocity tech¬ 
nique (3%) because while Rhaif changes with R 200 , the 
peak velocity occurs well within the radius at which the 
substructure was placed. 

3.4. Dependence on Particle Number 

Finally, we tested the performance of each technique 
for different resolutions (See Figure [8|. Beginning with 
the c = 10 halo containing 10 6 particles, transformed to 
be prolate with an e pro i a te = 0.3, random subsets of the 
particles were selected to create halos of the same size 
and shape, but lower resolution. 

Above 2 x 10 3 particles, all of the tested techniques 
are reasonably convergent (<10%) around the fiducial 
concentration for TV = 1 x 10 6 . For fewer particles, the 
residuals for the peak velocity technique climb quickly to 
> 80% of the fiducial value for 100 particles. TesseRACt 
and the half-mass non-parametric technique perform 
similarly on average, returning concentrations within 
20% of the fiducial value at the lowest resolution. How¬ 
ever, TesseRACt was more precise than the half-mass 
technique. Across the 10 halo realizations, the standard 
deviation in the concentrations returned by TesseRACt 
and the half-mass technique were 40% and 50% respec¬ 


tively at N = 100. Fitting was surprisingly accurate at 
low resolutions with residuals < 10% and standard devi¬ 
ations similar to TesseRACt. 

4. SUMMARY & DISCUSSION 

For idealized spherical halos, TesseRACt is slightly 
more accurate at recovering intermediate concentrations 
for N-body halos than techniques that assume spherical 
symmetry (^0.5% vs. ~1%). However, TesseRACt truly 
shines for non-spherical halos. For the most oblate or tri- 
axial galaxies, even the most accurate spherical technique 
using the peak circular velocity returned concentrations 
that underestimated the true value by up to ~10%, while 
TesseRACt had residuals of only ~0.5%. This is trou¬ 
bling given that halos in simulations are not often spher¬ 
ical. 

Studies of halos in cosmological simulations indicate 
that halos tend to be prolate overall, halos become in¬ 
creasingly triaxial at higher redshifts, and the most mas- 
sive halos at all redshifts are the least spherical (lAllgood 


_ ,_, . (Allg_ 

et al. 2006; B ett et al. 2007|fVera-Ciro et al.|2011| Despali 

et al. 2014p . Therefore, concentrations resulting from 


techniques that assume spherical symmetry would result 
in median concentrations that are biased overall to be 
lower than the true value and especially biased for halos 
that are massive or at high redshift. This bias could be 
expressed in the scaling relations between concentration 
and other intrinsic halo properties like mass, formation 
redshift, and environment. 

The concentration of a halo is b elieved to reflect th e 
state of the universe at its collapse ( Navarro et al.|1996 ) 


Therefore, on average, more massive halos should have 
lower concentrations, halos at higher redshift should have 
lower concentrations than those of the same mass at 
lower redshift, and halos in denser environments should 
have higher concentrations. However, the exact relation- 


vironment is still under debate dBullock et al. 2001 Eke 

et al. 2001 |Neto et al.' 2007 

; Gao et al. 2008; 

Klypin 

et al. 2011 Prada et al. 2012 

Bhattacharya et a 

l. 2013 

Dutton & Maccio 2014f Lud 

ow et al. 2014| Diemer & 

Kravtsov 2015). Since halo si 

rape also varies systemati- 
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cally with mass and redshift, it is possible that imposing 
the assumption of spherical symmetry could systemati¬ 
cally bias such studies. Use of a non-parametric tech¬ 
nique that is independent of halo shape like TesseRACt 
could be instrumental in eliminating this bias and pin¬ 
ning down these relations. 

In addition to biasing the median measured concentra¬ 
tion for halos in a given mass and redshift bin, the error 
introduced by assuming spherical symmetry would also 
increase scatter in the distribution of concentration for 
these halos due to scatter in halo shape. Intrinsic scatter 
in this relationship is expected to result from differences 
between any two halos’ mass assembly histories that re¬ 
sult in slightly different concentrations. However, if scat¬ 
ter is also being introduced by the technique used to mea¬ 
sure concentration, measurement of this scatter becomes 
less informative. Since scatter in the concentration-mass 
relationship is greatest for those halos which are t he least 


virialized and ha ve the poorest NFW fits (Jing 2000 


Neto et al. 2007), removing these halos can help to re¬ 


move some of this bias. However, even virialized halos 
can be non-spherical and it is likely that the assumption 
of spherical symmetry would still contribute to the re¬ 
maining scatter. TesseRACt could help to identify what 
the true intrinsic scatter in the concentration-mass rela¬ 
tionship is for both virialized and un-virialized halos. 

TesseRACt could be improved for future studies in a 
number of ways. One possibility is using tessellation in 
the full six dimensional phase space to place better con¬ 
straints on how deep a particle resides within the host 
galaxy’s potential. This would be particularly useful for 
correcting for substructure. While substructure particles 
may be confused with parent halo particles on the basis 
of spatial volumes alone, the two would be more easily 
differentiated in velocity space. It may also be possible 


to improve TesseRACt’s performance by more accurately 
parameterizing the relationship between traditional con¬ 
centration and that calculated using tessellation based 
volumes as in Eqn. [8] Fitting over a wider range of 
concentrations and allowing a more complex relation¬ 
ship should improve the slight dependence of TesseR¬ 
ACt’s performance on concentration. However, since 
this dependence only affects concentrations by < 3%, 
the parameterization presented here should be sufficient 
for most studies. 

Although only spatial information was used here, tes¬ 
sellation can also be used to glean information non- 
parametrically from simulations by tessellating over 
other parameters. In simulations which include gas, tes¬ 
sellation over parameters like star formation rate, tem¬ 
perature, and metallicity in addition to position and ve¬ 
locity can help to identify structure and explore forma¬ 
tion scenarios. 


5. CODE 

The TesseRACt method for estimating concentrations 
from Voronoi volumes has been assembled into a publicly 
available Python package that will be made available at 
the time of publication. This release will include example 
halos, routines for running a nd handling output fro m an 
adapted version of VOBOZ (|Neyrinck et al.l|2004|), and 
routines for estimating concentration using the spherical 
techniques tested here. 

K.H-B. acknowledges the support of a National Science 
Foundation Career Grant AST-0847696 as well as the 
supercomputing support of Vanderbilt’s Advanced Cen¬ 
ter for Computation Research and Education. Support 
for M.L. was provided by a National Science Foundation 
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